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Abstract 

The Stieltjes constants 7fc(a) are the expansion coefficients in the Laurent 
series for the Hurwitz zeta function about its only pole at s = 1. We present 
the relation of 7fc(l) to the rjj coefficients that appear in the Laurent expansion 
of the logarithmic derivative of the Riemann zeta function about its pole at 
s = 1. We obtain novel integral representations of the Stieltjes constants and 
new decompositions such as S^fn) = Sy(n) + S\(n) for the crucial oscillatory 
subsum of the Li criterion for the Riemann hypothesis. The sum S^(n) is 
0(n) and we present various integral representations for it. We present novel 
series representations of ^(n). We additionally present a rapidly convergent 
expression for 7^ = 7fc(l) and a variety of results pertinent to a parameterized 
representation of the Riemann and Hurwitz zeta functions. 
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Introduction and statement of results 



In the Laurent expansion of the Riemann zeta function about s — 1, 

\ n 



i °° c_iv 

coo = — 7 + E - ^ 

s — 1 ~i n! 



n=0 

the Stieltjes constants 7& [91 [TO], [22j (25J [271 E2] can be written in the form 



x 1 k \n k+1 N' 



7 fc = lim ( V — ln fc m - — — ) . ( 2 ) 

1 iv-,oc^ m k + 1 1 



From the expansion around s = 1 of the logarithmic derivative of the zeta function, 

COO i 



oo 



C( S J - - P =o 

we have 



l-Y,Vp(s-V P , |s-l|<3, (3) 



lnC(s) = -ln(s-l)-^^(s-l) p . (4) 

p=i P 

The constants 77^- can be written as 

(-1)* / " I . ln fc+1 iV\ 
% = ,, hm > — A(m) ln fc m — , (5) 

where A is the von Mangoldt function [T5], [22J [2U, [33l [30], defined by A(n) = hip if 
n = p h for a prime number p and some integer > 1, and A(n) = otherwise. The 
radius of convergence of the expansion (3) is 3, as the first singularity encountered is 
the trivial zero of ((s) at s — —2. The Dirichlet series corresponding to Eqs. (3) and 
(4) valid for Re s > 1 are 

CM = _£ *2> and k fl,] = £M. (6 ) 
Indeed, Eqs. (6) hold on a = Re s = 1 if t = Im s 7^ [33] . 



We note that there are various explicit expressions for the von Mangoldt function. 
These include 

LCM(l,...,n) 

6Xp[A(n)] = LCM(l,...,n-l) ' (7) 
where LCM is the least common multiple function, and Linnik's identity [28J (pp. 

21-22) 

A(n) Wln2/_ i)k 

In regard to Eq. (8), the strict r' k and exact r k divisor functions are related by 

T k( n ) = \{ n ii ■■■,n k >2;n 1 ---n k = n}\ (9a) 



Wt(n). (96) 



In particular, there is a finite sum on the right side of Eq. (8). 

We recall that the function In ((s) is intimately connected with the prime counting 
function 7r(x), the number of primes less than x. We have 

\n((s) = s T -?P—dx. (10) 

J2 X{X S — 1) 

Hence the behaviour of the function n(x) is related to the important coefficients 
r]j. For further background on the classical zeta function we refer to standard texts 

H51I221I2SG31I3D1. 

The Hurwitz zeta function, defined by £(s, a) = I]^L (n + a)~ s for Re s > 1 and 
Re a > extends to a meromorphic function in the entire complex s-plane. The 
generalization of the Laurent expansion (1) is 

1 00 ( — 1 \n 

C{s, a) = - + E ^7n(a)(, - 1)". (11) 



As shown by Eq. (1), by convention one takes 7^ = 7^(1). 

Having set the notation above we may now state our main results. 
Proposition 1. For integers k > we have 

(-l) fc (-1)" ^ [1-A(m)l k , 10 . 

lk-r] k = ^-rj- J2 ' ln m ' ( 12 ) 



k\ k\ ^ m 

where the sum on the right starts with m = 2 when k > 0. The case of k = yields 
the identity 
Corollary 1 

7= lf;[lz^01, (13 ) 

Z m=l ' Ai 

where 7 is the Euler constant. 
Put 

S,{n) = -£ ("Wi- (14) 

m=l V'V 

Let L%(x) be the Laguerre polynomial of degreee n (e.g., [HE]) and -Pi(t) = -&l(£ — 
[t]) — t — [t] — 1/2 the first periodized Bernoulli polynomial (e.g., [22], [33] ) . Then we 
have the series representation given in 
Proposition 2. For integers n > 1 we have (a) 

S 2 (n) = S 1 (n) + S A (n), (15) 



where 



S» = t ^ ftVa = f - t L MdPl (t), (16) 



™ [1-A(m)] rl 

m=l 771 



n + S 2A (n) = n + £ [1 A(ro)1 £^(hm), (17) 

^2 m 



(b) S^n) = O(n), (c) the average values 



1 1 r 00 1 

mE S t(») = m I 1 Ll_ 1 (\nt)dP l (t), (18a) 
T, £ «*(») = S ^ I ^ M ^-.0™). < 186 ) 

JW 71=1 iW m=l '' t 

and (d) for iV > 1 a fixed integer 

S» = - E ^-^^ - L n (ln TV) + 1 + ^L^lnN) + (^_) 

where e > is arbitrary. 

We have found new integral representations of the Stieltjes constants, given in 
Proposition 3. Let Re a > 0. Then we have (a) 

/ \ 1,* ln fc+1 a , 2 /-°° (y/a - i) ln fc (q - iy) 
7fc(fl) = Ya ln a ~ TTT + a Re Jo (l +y yai)(e^-l) dy - (19) 

(b) Let n > 1 be an integer. Then we have 

_ ^ ln fc (m + q) ln fc (n + a) _ ln fc+1 (n + a) 

2 y°° [y/(n + a) - i] ln fc (n + a- jy) 

n + a io [l + y 2 /(n + a)*](e**v-l) V ' 1 J 



(c) 



7r X roc \ n m+1 (2a-l-it) , 
7m(a) = -77 7T Re / 1 27 77777 rft 



+![ m ! V "V, v v -Re / ,\ " v, dt. (21) 

2 h j\ (m-j + l)\ Jo cosh 2 (7rt/2) V ; 



m+l 



2 m + 1 Jo coslr(7rt/2) 
ln j 2 (-iy +1 /•~ln m - J+1 (2a-l-^) 
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Let if>(z) = T'(z) /T(z) be the digamma function, where V is the Gamma function 
[H HI EE]- From Proposition 3 follows Corollaries 2 and 3. From the case m = in 
Proposition 3 we have 
Corollary 2. (a) 

"I 2 r°° (via i) 

7oW = -^(a) = --l n0+ -Re j ( ' A,, (22) 



(b) for integers n > 1, 



n ^ 1 1 , / 2 r. f°° \y/(n + a)-i] 

7o (a) = -^(a) = E^+^^y-ln^ + ^^yRe ^ p + ^ - _ ^ dy, 

(23) 

and (c) 

, / x r°° \n(2a - 1 - it) , 

4 Jo cosh 2 (7rt/2) V ; 

Based upon the a = 1 case of Proposition 3 we obtain additional new integral 

representations of the sum 6" 7 defined in Eq. (16): 

Corollary 3. 

f°° 1 r i 

(25) 

and 

S» = - / [L n (2a - 1 - it) + L n (2a - 1 + it) - 2} 

4 jo cosh (irt/2) 

7r " ln j 2 /n\ z 100 (it 
- 7 E— r ■ / {i^iL7'-n;j>l;ln(2a-l-it)]+ 1 F 1 [ 7 --n;j+l;ln(2a-l+it)]}— -j- 
4 ~[ j! \J/ Jo cosh (7 



cosh 2 (7rt/2) 
(26) 



where 1-F1 is the confluent hypergeometric function, with (^j \Fi(j — n;j + l;x) = 



L J n -j(x). 



Let d(n) be the number of divisors of n. Then we have 
Proposition 4. (a) We have for integers m > 
00 ln m n r 



^ — [d(n) - Inn - 27] = (-l) m+1 ( 1 + -^7 ) 7m +i - 2 7 (-l) m 7r i 
~[ n V m + 1 / 



+E(-i)'(7j 7 ^- (27) 

In particular we have at m = 
Corollary 4. 

E -[d(n) - Inn - 2 7 ] = -3 7l - 7 2 = -3 7l -iff H^MV (28) 
giving 

Corollary 4'. 



~ d( m ) n , 

> — — - - hm 

m=l 



In n + 2 7 In n 



7 2 - 2 7l > 0. 



(b) We have for integers m > 



00 i n m 2( — l) m+1 

v"l) m E - A M Inn - 2 7 ] = (m + l) Wl + 1 . 7m+ i - 2 7 (-l) m 7 „ 

~[ n m + 1 



+ E(-i)^;j7n^- (29) 

In particular we have at m = 
Corollary 5. 

00 j 

- [d(n) - A(n) In n - 2 7 ] = 771 - 2 7i - 7 2 = 0. (30) 

71=1 ^ 



We may obtain new series representations of S y (n) from part (a) and of S^n) from 
part (b). As an illustration we have 
Corollary 6. For integers n > 1 we have 

S 2 (n) = 7 n+E Li ^ nm \ d{m)-K{m)\nm-2 1 }-2±{-ir( n \- ^ 

Let n{n) be the Mobius function. 
Proposition 5. Then we have (i) 

Vo = -l = l±^^k, (32) 

Z fc=2 K 

and for integers £ > 1 

<-^-^£^* + £^£^» (33, 
(ii) For integers j > 1 we have 

* + e + i! 7 ' f h n\ 7 "(i-»)! ln 



-n+1 



From Proposition 5(i) follows 
Corollary 7. We have 

*<») - E^m, [^(m*) - ij + E — B-i) y g^r^ 

(34) 

Proposition 6. Let £?, represent the Bernoulli numbers. For integers m > 1 and 
Re A > we have 

m! A 5 m _ m ln m " £ 2^/ A \ fe * n , j fk\ 1 h/(j + 1) 



7m 1 + A§(m -€ + 1)! i\ £(l + Aj g ( ^ 1^ (j + 1) 



k k /i„\ l„w+l, 



i i a Vvv ^,Y^in m+1 (j + i) fi -+i , lnm+ i 2 . 



(l+A)ln2(m + l) fe fA 1 +^/ jti \j) O' + l) 

(35) 

For m = we have 

1 1 ^/ A \ fc * N ,/A;\ln(j + 1) , . 

Proposition 7. (parameterized series representation of the Hurwitz zeta func- 
tion) For Re s > 1, Re a > 0, and Re A > we have the representation 

c(s ' 0) = |(r^gC)^(JT7?- (37) 

Proposition 8. Let the polylogarithm function Li s (z) = J2T=i z k /k s = z<&(z, s, 1) 
where $ is the Lerch zeta function, s G C when \z\ < 1 and Re s > 1 when \z\ = 1. 
Then for g > 1, Re t > 0, and m > an integer we have the integral 

I ^ - *>'"" Li << 1 - - „t (^TTF^- 1 ^ Q ln(t + A (38) 

giving 

Corollary 8. At q — 1, with Lii(^) = — ln(l — 2) we have the special case 

-B(^)=tj^^t(-^)H^), (39) 

where B(x,y) = T(x)T(y)/T(x + y) is the Beta function [I]. 
Let for Re x > and Re y > 

B( Zl ,z 2 ,x,y)= re- l {l-ty- l dt (40) 

J Zl 
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be the generalized incomplete Beta function and Ei(x) = j^ 00 {e t /t)dt the exponential 
integral. Then we have 

Proposition 9 (context of the logarithmic integral) 



— (1 ~ «r 'du = ^2(-iy [ V . ){Ei[(a + j)\nz l ]-Ei[(t + j)lnz 1 ] 
J xi mw j_Q \ 3 J 

-Ei[(a + j) In z 2 ] + Ei[(t + j) In z 2 ]} , (41) 



giving 

Corollary 9 

/ -du = Eifalnzi] -Eiftlnzi] - Ei[aln^ 2 ] + Ei[t) \nz 2 } 

J Zl In it 

-Ei[(a + l)lnzi] + Ei[(* + l)lnzi] + Ei[(a + 1) lnz 2 ] - Ei[(i + 1) lnz 2 ]. (42) 

Propositions 1 and 2 were first obtained by the author several years ago as com- 
plements to Refs. [8] and [TT] concerning the Li criterion for the Riemann hypothesis. 
The proof of Proposition 2 given in the sequel is more direct than the original. The 
sort of alternating binomial sums that occur in Proposition 6 motivates a study of in- 
tegrals such as appear in Propositions 8 and 9. After the proofs of these Propositions 
we present Discussion that contains additional examples and extensions. 

Proofs of the Propositions 

Proposition 1. In order to derive Eq. (12), we add Eqs. (1) and (3) and make use 
of Eq. (6), resulting in 



-if 

In - Tin 



nl 



;._ ir = f;M, (43 , 

k=l K 
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wherein the pole at s = 1 has been removed. Repeated differentiation of Eq. (43) 
gives 

Corollary 10 



C W 00 + 



COO 

C(s) 



(J) 



E 

ra=j 



:-i) n 



n! 



"7n - «n 



n(n - l)(n - 2) • • • (n - j + l)(s - i; 



n— j 



^iFE' 1 -.^ 1 .^/.. 



fe=l 



(44) 



where the sum on the right starts with k = 2 when j > 0. Taking s — > 1 + in Eq. (44) 
gives Eq. (12). 

Remarks, (i) Taking different values of s in Eq. (44) yields various connections 
between the 77's and the Stieltjes constants. For instance, with s = 2 in Eq. (44) we 
have 



E 

n=j L 



4)n 



-7, 



n ' In 



-n) 



k=i 



k 2 



(45) 



where (a)p = T(a + i)/Y(a) is the Pochhammer symbol, (ii) For Eq. (13) we use the 
values 70 = 7 and 770 = ~ 7 (e.g., [8], [Tl] ) . This special case (13) for the Euler constant 
has unfortunately appeared in several places in the literature including [16] and [19] 
(p. 109) with the sum starting at 2 instead of 1, missing a dominant contribution of 
1/2. (iii) Proposition 1 may also be proved by forming 

[AM - 1] 



In COO 



[CO/) - l]dy = £ 



7 - c w + £ 



where [5] (p. 156) 



C\ 1 = lim 



E 

.n=2 



^-1 



1 

ralnra 



+ 



n=2 

4)i+i 



n s Inn 



7i-i 



1) J 



(46) 



In (In 2;) 



0.7946786. 
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The one new sum that appears here has value 

f^")- 1 ]^-^ -0.217464. 

(iv) The Stieltjes and rjj constants are strongly connected and one may easily write 
a recursion relation between the two sets of constants [8 J (Appendix). For instance, 
rji = 7 2 + 271. Then from Eq. (12) we may write either 

1 2 2 ~ [1 - A(m)] . 
3 m 



or 

1 

7i = o 



~ [1 - A(m)] . 2 



with the corresponding formula for 7 given in Eq. (13). Similarly, we have corollary 
expressions for all of the 7^ and rjj constants. 

Proposition 2. By applying Proposition 1, the definition (14) of ^(n), and the 
power series definition of L™, we obtain 

*(») = t(^Ji(3^. + S*(»). (48) 

The second line of Eq. (17) follows since A(l) = and L^^O) = n. We obtain 
alternative forms of the sum S 1 {n) of Eq. (16) by using the integral representation 



/oo \ poo ]ti 1 + 

-ln k tdP 1 (t) = ^ — ^— (k-ln^P^dt- 8 k0 /2 (49) 



7fc 



■D^- 1 /■<» /c/\ fe ln fc t 



(yfc - 1)! Ji RW U* 



p fc(*)hl ^dt-5 k0 /2, (50) 
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where 6jk is the Kronecker symbol. In Eq. (50), where we integrated by parts k — 1 
times, Pk(t) = Bk(t — [£]) with Bk the /cth degree Bernoulli polynomial, such that 
P' k+ i(t) — {k + l)Pfc(t). By then applying the definition (16) and the power series 
form of the associated Laguerre polynomials we have 

Pi(t) 



t 2 



L 2 n _ 2 (\nt) + L^ilnt) dt + n/2. 



(51) 



This representation of S y is equivalent to integration by parts on the expression given 
on the right side of Eq. (16). 

For part (b) we observe P\{t) =0(1). So for a constant C > the integral term 
in Eq. (51) is majorized by 



L 2 n _ 2 (lnt) + L^lnt) 



C 1 + 



dt = C I e 
'o 

1 



n-2, 

where we used the value of a Laplace transform [T7j 

1 (n-l)\ 



Ll_ 2 (u) + Ll_ 1 (u) 
n > 3, 



du 



e~ x U 



lie j djiXs 



(52) 



(53) 



/o T(z/) (n — v)\ 

Then S 1 {n) = 0(n). For part (c) we apply Eqs. (16) and (17) and interchange the 
order of operations. 

For part (d), we make use of [27] (Section 2) 



^ In* v ln fc+1 n 1 In* n 



V 



jfe + 1 2 71 



+ 



1 



2-e 



(54) 



with e > arbitrary. This equation is essentially an asymptotic form of 7& for k 3> 1. 
We substitute this equation into the definition (16) of S 7 (n), apply the power series 
form of the associated Laguerre polynomials, and Proposition 2 is completed. 
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Proposition 3. Part (a) is based upon the well known Hermite formula for £(s, a): 

10 sin(s tan -1 y/a) dy 



a a 



(y 2 + a 2 ) s / 2 (e 2n y-l)' 



We use the expression 



! 1 

tan u> = — In 

2i \ 1 — iw 



(55) 



(56) 



so that exp[±istan 1 (y/a)] = [(1 + iy/a)/(l — iy/a)] ±s ^ 2 . Then the integral term in 
Eq. (55) becomes 

sin(stan _1 y/a) dy 



7 

Jo 



(y 2 + a 2 ) s l 2 (e 2 ^ - 1) 



Jo 



1 



(y 2 + a 2 ) s / 2 
1 



1 + iyja 
l-iy/a / 



s/2 



1 — iyja 
1 + iy/a, 



s/2 



dy 



a Jo 



a s Jo {l+y 2 /a 2 ) 

\y/a-i)e 



(1 + iy/a) (1-iy/a) 



(1-iy/a) 3 - 1 (1 + iy/a) 3 - 1 



(e 2 *y - 1) 
dy 

[e 2 *y - 1) 



(l + y 2 /a 2 ) 



-(s-l) \n(a-iy) 



+ (y/a + i)e 



-(s— 1) ln(a+iy) 



dy 



(e 2 ^ - 1)' 
(57) 



For the other terms in Eq. (55) we have 



2 + s-l 



q1— s g— (s— l)lno g— (s— l)lna 
+ 



2a 



1 1 00 f-lV 



S-l 
21 (-l) j+1 



ln J+1 a(s - 1) J 



(58) 



j=0 J- j= 

We then expand the right side of Eq. (57) in powers of s — 1, combine the result with 
Eq. (58), and compare with the defining expansion (11) for 7fc(a). We find 
1 



lk{a) = —\n K a- 



k + 1 a Jo 



(l+y 2 /a 2 )(e 2 ^- 1) 



(59) 
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that is the same as Eq. (19) and part (a) is proved. 

For part (b) we apply Abel-Plana summation (e.g., [31 J p. 90) to write for all 
complex s ^ 1 



n-l 



^ 1 (n + a) s (n + a) L 



+2 



°°sin[stan 1 y/{n + a)} dy 
o [y 2 + (n + a) 2 ] s / 2 (e 2 ^-l)' 

(60) 



We then expand the terms of this equation in powers of s — 1 in like manner to part 
(a). 

Part (c) makes use of the integral representation for Re a > 1/2 [23] 



CO, a) 



[8-1) 



[t2 + (2a _ 1) 2 ] (l- S )/2 



cos[(s - l)tan~ 1 [t/(2a - 1)] 



cosh 2 (vrt/2) 



dt. (61) 



We again apply Eq. (56) so that we are able to write 



CO, a) 



7T 1 



1 ^ ln ' 2 / IV f 



4s 



+ 



(2a - 1 - it) 8 - 1 (2a - 1 + itf- 1 



7T 1 

4~ 



- / y h/(2a - 1 - it) + h/(2a -1+it) (s- 1)' 



dt 

cosh 2 (vrt/2) 
dt 



cosh 2 (vrt/2) 



7T 1 ~ In 3 ' 2 ~ (-1)' /-oo r 



4s 



x 3 =i J- e=o 



t\ Jo 



h/(2a - 1 - it) + ln € (2a -1 + it) (s-1) 



cosh 2 (7rt/2) 



(62) 



We then (i) separate out the £ = term two lines above and use the integral 
Jq°° sech 2 (7rt/2)dt = 2/ir, and (ii) reorder the double sum in the last line of Eq. (57), 
thereby obtaining 



CO, a) 



1 7T 



-+-Y 

- 1 - ^ e=o 



'1) 



+ 1)! 



ln m (2a - 1 - it) + h/ +i (2a -1+it) 



e+i, 



dt 

cosh 2 (vrt/2)' 
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oo m+j 

+7E(-i) m E 



m=0 j'=l 




\n j 2 f 
j\ (m — j + 1)! io 



\n m - j+1 (2a - I - it) + \n m ~ j 



(2a - 1 + it) 



(63) 



x 



cosh 2 (vrt/2) 



Comparing Eq. (63) with the Laurent expansion (11) gives part (b) and Proposition 
3 is proved. 

For Corollary 2 we put m = in Eqs. (19), (20), and (21) and use the fact that 



Remarks. As it must, Eq. (60) satisfies ^(0, a) = 1/2 — a, ('(0,a) = lnT(a) — 
(1/2) ln(27r), and B m (x) = —m((l — in, x) for positive integers m. 

Equation (22) recovers the result of differentiating Binet's second expression for 
lnr(2i). Equation (24) subsumes the special case at a = 1/2 given in [T7] (p. 580). 

For Corollary 3 we apply the definition (16) of 5 7 (n) together with 
Lemma 1. We have 



that follows from the power series form of the associated Laguerre polynomials. The 
second line of Eq. (26) follows first from reordering a double sum to obtain 



7o(z) 



■1>(x) (e.g., 0). 




(64) 





(65) 
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We then use various relations to express in terms of Pochhammer symbols: 



n 



k + j (k+j)\ 



-l) k +i (-l) k+j (-l) k+j C-lVn' 



(fc + j)! 



(&+j)! (n-i)! 



(66) 



where + j)! = + 1)*. We then apply the power series definition of the confluent 
hypergeometric function, Eq. (26) follows and Corollary 3 is completed. 

Remarks. Similar to Proposition 3 we have obtained further integral representa- 
tions of the Stieltjes constants and the sum S 1 (n) from related integral representations 
of the Riemann zeta function. In particular, we have [23] 

2 poo (J^t/ 



as) 



--2 s (l + t 2 )- s / 2 sin(stan- 1 t)^— . 

i Jo v ; v V* + r 



and 



As a byproduct we have from Eq. (68) 
Corollary 11. For all complex s ^ 1 we have 



_ s / 2 cos(stan 1 1) 



cosh(7rt/2) 



dt. 



(67) 



(68) 



C(s) 



2(1 - 2 



l-s\ 



+ 



1 



1 



{1/2 -iw) s (1/2 + iw) 1 
1 dw 



dw 



cosh(7ru>) 



(69) 



2(1 - 2 1 - 3 ) J-oo (1/2 - iw) s cosh(vrw) ' 
Equation (69) follows from Eq. (68) with a simple change of variable and the use of 

relation (56). 

Of a slightly different flavor, we have used the representation [Ej valid for < 
Re s < 2 



CM 



1 



sin us 



S — 1 7t(s — 1) JO 



11/ (t) 



1 + t 



(70) 
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to obtain 

Corollary 12. We have for integers m > 



[m/2] 



-l) k+1 n 2k 



'-r) m m \ v 

^ ) ■ ^ (2k + l)\(m-2k)\Jo 



In 



m—2k 



11/ (t) - 



1+t 



dt. 



(71) 



The proof of Corollary 12 makes use of the Taylor series 



sin tcs 



E 



:7T 



2fe+l, 



8-1) 



2fe+l 



(72) 



fc=0 (2* + l)!' 

together with series manipulations. The case of m = in Eq. (71) recovers 70 = 7, 
as is readily seen by a limiting argument. Corollary 12 provides another form of the 
sum Sj(n), that we omit. 

We have found that Proposition 3(a) subsumes the a = 1 case derived in Ref. 
[2] by a contour integration. This reference evidences that Proposition 3 should be 
a practical method for computing 7 fc(a). Indeed, there is now an arbitrary precision 
Python implementation [29] . 

Proposition 4- For part (a) we form the combination of Dirichlet series 



( 2 (s) + C(s) - 2 7 C(s) = E ~[d{n) -Inn- 2 7 ] 

~i n 



(73) 



n=l 



wherein simple and double polar terms have been eliminated. By the use of Eq. (1) 

and series manipulations we find 

00 ^ 

J2 -[d( n )-lnn-2 7 ] 



n=l 
00 

= E 

j=0 



rr 



2(-iv +i (-iy +1 i-iy ^ i-iy 



itij-i 



(s-iy. 

(74) 



Taking the limit as s — > 1 + in this equation gives Corollary 4. Corollary 4' then 
follows from the limit relations (2). Taking m derivatives with respect to s in Eq. 
(74) and putting s — > 1+ yields 
00 \n m n 

(-l) m 5] [d(n)-lnn-2 7 ] 



n=l 



n 



i \m+l / 0\ ('_1 N \ m m 

' 1 + ^ 7m+l - 2 7 ^-^7m + E 



m! 



£!(m-£)! 



lllm-t 



ml '•■ m+1 
and this gives Eq. (27). 

For part (b) we form the combination of Dirichlet series 

(\s) - (y) (s) - 2 7 C( S ) = £ ~ s [d(n) - AH Inn - 2 7 ]. 



m\, (75) 



c 



71=1 



(76) 



We then use both Eqs. (1) and (3) and series manipulations to find 

CO J 

E — [d(n) - A(n) Inn - 2 7 ] 



n=l 

oo 

= E 

J"=0 



^ ■ 7j+1 + (j + - 2 7 ^-^ 7j + E ^, 7£7,-£ 



(J + 1)! 



{s-iy, |s-i|<3. 

(77) 



Taking the limit as s -> 1 + in this equation gives Corollary 5. More generally, by 
taking m derivatives of relation (77) we have 



00 ln m n 

-l) m E^[ d W- A (™) ln ™- 2 7] 



n=l 



n" 



= E 



2{-iy +1 



-1)' 



yr^+i + (J + l)^+i " 2 7 L -y^ 7 , + E 



[-1Y 



xj(j-l)(j-2)--.(j-m+l)(s-iy- m 



lelj-e 



\s-l\ < 3, 
(78) 
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wherein term-by-term differentiation is justified within the stated radius of conver- 
gence. Taking s — ■> 1 + in this equation gives Eq. (29). 

For Corollary 6 we write 77m— a from Eq. (29) and apply the definition of 5*2 of Eq. 
(14) in the form S^n) = ^n — J2m=2 (m) 7 ?™- 1- ^ e then apply Lemma 1 for the power 
series form of L l n _ x to obtain Eq. (31). 

Proposition 5. (i) We begin by writing 



COO 

as) 



-E 



A(ra) 



Re s > 1, 



n=l k=l 

and forming ('(s) from Eq. (1). This gives 

COO ^ M*) ^ (- 1 )' 



COO 



y-^ M^O 

fc=2 ^ £=1 



■ln'JfeOj-l) 



L ) j=0 J- 



(79) 



(80) 



where we have used J2kLi n{k)/k = 0. We then multiply the series, reorder a double 
summation, and make use of 



fc=2 



A- 



1. 



wherein we appealed to a Tauberian theorem ([34J, Ch. V). After these operations we 
compare with the defining Laurent expansion (3) for the rjj coefficients, and Proposi- 
tion 5(i) follows. 

For part (ii) we instead use the identity 

A 1 

(82) 



COO ^" J ds(( s y 

We use the Laurent expansion (3) for the left side, and expansion (1) for £(s) and 
the Dirichlet series for 1/C(s) on the right side. Expanding in powers of s — 1 gives 



1 3=0 



1) J 



7 + £ 

-L ^ — n 



■l) n 



n=0 



7n(*-l) n 
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00 C— lV 

1 + 53 ^p- ln ' ^ - !) J 

i=i J - 



J3) 



We then expand the right side, reordering the last series, apply relation (81) and (ii) 
follows. 

Corollary 7 follows from part (i) by using the definition (14) of S^n) and applying 
Lemma 1. Similarly, another form of $2(71) could be written based upon the expression 
for rjj given in part (ii). 

Proposition 6. The key starting point of the proof is the Amore representation of 
the Riemann zeta function [3] 

(84) 

Beyond the condition Re s > stated in Ref. [3], the representation (84) holds for 
all complex s 7^ 1, as the summation continues to converge for Re s < 0. (We further 
discuss Eq. (84) in the following section.) So Eq. (84) is globally convergent, as is its 
A = 1 special case embodied in the Hasse representation [18], [9] . Moreover, beyond 
the original condition A > of Ref. [3], we may take A complex with Re A > 0. 

The proof now proceeds as in Proposition 6.1 of Ref. [10J. However, that descrip- 
tion was terse and we now have the (arbitrary) complex parameter A. Therefore, we 
believe it is worthwhile to supply some more details. We first write again 

U + ir = (j + l^expf-lnO" + - 1)] 

1 °° c_iy 

= -t^j: [ -4 l ^(j + i)(s-iy. (85) 
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From the generating function of the Bernoulli numbers, we have 

f„t j- oo in oo rn 

' =E^(1)^ = E(-1)^A, \t\<2n, 



e* - 1 1 - e 



n=0 n - n=0 

where B n (x) is the nth Bernoulli polynomial. We then put t — (s — 1) In 2 in Eq. (86) 
and obtain 



(86) 



(1-2 



l-SN-l 



|.-1|<£, (87, 



ln2( s -l) ^ + 1)1 
We substitute Eqs. (85) and (87) into (84), writing the sum over t in Eq. (84) as 
the I = term and the rest of the terms. We use 



k k 



: \ x + ^/ j=0 



E(-l) J nU7- 1 TT = (A + l)ln2, 



obtaining 



C(*) 



1 + A 



vJ7 * U + 1) 



(88) 



In2( a -1) § (j + 1)! 



hy'2( s -i) j 



x 



fc=o \ x ^ 'V j=o 



59) 



We multiply the terms in Eq. (89), separate out the simple polar part, and compare 
with the defining expansion (1) for the Stieltjes constants. For the last product of 
series in Eq. (89) we use the reordering 

oo oo oo oo oo m 

E EC ■■)(«- i) n+ ' = E E (■■■)(« - !) m = EEC- •)(« - i) m , (90) 

n=0£=l l=\m=l m=ie=l 



and the Proposition follows. 
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In particular, at A = 1/2 we obtain 



Corollary 13. 




ni 



(m-e + iy. t\ 



Bm-e+i In' 




31n2 (m+ 1) 



2 1 




k\ ln m+1 (j + l) 
3 J U + 1) 



ln m+1 (A;+ 1) 



(m + 1) 



In 



m+l 



2. 



(91) 



Remarks. The expression (91) may be attractive for some computational appli- 
cations because it exhibits even faster convergence than the form of the Stieltjes 
constants resulting from the Hasse representation of the zeta function |10j . The free 
parameter A in Eq. (84) may even be taken as a function of s. Based upon the 
principle of minimal sensitivity, one chooses to lowest order A 1 - 1 ** = 2 _s . Amore [3] has 
noted that with this choice in Eq. (84) a still exact expression for the zeta function is 
obtained. Since in obtaining the Stieltjes constants we expand about s = 1 we chose 
A = 1/2 in Corollary 13. The form of ((s) given in Eq. (23) of Ref. [3], where s 
dependence now appears in 4 places, could again be expanded about s = 1 but this 
additional effort does not seem warranted. 

Proposition 6 shows once again that the Stieltjes constants may be written in terms 
of the Bernoulli numbers and elementary constants such as powers of the logarithm 
of the natural numbers. In turn, the same holds for any quantity expressible in terms 
of the Stieltjes constants. While not many arithmetic properties are known for the 
Stieltjes constants, the Bernoulli numbers have many known important properties 
such as congruence or other relations. For instance, writing the (rational) Bernoulli 
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number Bj = numer(Bj)/denom(Bj), we have denom(i? 2r) ,) = ri(p-i)|(2n)P where the 
product is over prime numbers p such that p — 1 divides 2n. 

Proposition 7. The Hurwitz zeta function has for Re s > 1 and Re a > the 
integral representation 

C(s,a) = / dx = — — / v 7 / J — tte. 92 

sv ' ; T{s)Jo e x -l T{s)Jo (l-x) K J 

Then we introduce a parameter A with Re A > 0, make use of a geometric series 

expansion for x G [0, 1], and follow this with a binomial expansion: 

1 



C(s,a) 



x a - x 


"^(l/x)] 5 - 1 


(1 + A) 


>-(m)l 



1 oo \fc fc 1 -i 

= J_V— A— rE ^ / x^'-Xl/a;)]- 1 ^. 
r(a) fro (! + X ) k+1 U W Xl Jo 

Performing the integral in terms of the T function gives the Proposition. 

Proposition 8. We first obtain an alternating binomial sum by integrating the 
Beta function. We have 

B{x,y) = f 1 u x - 1 {l-u) y - 1 du = 2 r /2 sin 2:E - 1 0cos 22/ - 1 ty = B(y,x), 
Jo Jo 

min[Re x, Re y\ > 0, (94) 



k 



dx 



(93) 



so that 



f B(x,y)dx= f 1 ^ \ Ua 1) {\-uy- 1 du (95a) 
J a Jo mw 



3=0 
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jr(-iy( y . ^Kt + jO-Ma + j)], (95c) 



3=0 



J 



where the form (95b) follows by binomial expansion in Eq. (94). Upon comparing 
Eq. (95a) and (95c) we have 



(l- u y- 1 du = J2(-iy( y . m(t + j)+C, (96) 
j=o V 3 J 



o In u 



where C is a constant to be determined. A simple way to do this is to put y = 2 
whereupon Jq 1 -u t_1 [(l — u)/(\nu)]du = lnt — ln(l + t). Therefore, C = and 



[ 1 ^(i- u y- 1 du = f:(-iy( y . ^Ht + j) 

Jo In « p Q \ 3 J 



(97) 



a result closely related to tabulated integrals when y is an integer [T7] (p. 546). Now 
if y — n + 1, n > 1 an integer, the sum in Eq. (97) terminates: 

r T— (1 " = E(- 1 ) J ("1 + j)- (98) 



3=0 



We next multiply each side of this equation by l/(n — in + l) q and sum on n from m 
to oo. We shift the summation index on the left side, apply the series definition of 
the polylogarithm function, and the Proposition follows. 

Remark. The Proposition may be extended by analytically continuing to Re q > 1. 

Proposition 9. By applying binomial expansion to Eq. (40) we have 



B(zi, z 2 , x, y) = E(-l) j P ^7 J . (99) 



3=0 



We then integrate this expression on x from a to t and compare with the integrated 
form of Eq. (40), 

B(z 1 ,z 2 ,x,y)dx = - -(1 - uY^du, (100) 

Jzi In u 
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and the Proposition follows. Corollary 9 is the special case of y — 1. As a further 
special case we have 
Corollary 14 

rl rz 2 (I — U a - l \ 

/ B(2,z 2 ,x,l)dx= / ^— '-du, (101) 

J a J 2 In U 

showing the close connection with the logarithmic integral Li(z) = f£ dt/lnt. 

Discussion 

Amplifying that the representation (84) holds for all complex s ^ 1 we easily 
verify that 

1 00 / X _ 1 \ fc 1 

w-tta S(xn) - 5 - « 

Furthermore, we have 



C'(s) = - 



1-s 



1 A \ 1 ' mi ln(l+j) 



(i + a)^ vi + a; ^ v^/ AJ u + 1 ) 



s G C, s ^ 1, (103) 



so that 



A \ fc * ,,AA 1 



^"■"^(hooSIttaJ g ( - 1)4 Uv k(1+ * (104) 

This value is easily shown to be —(1/2) ln27r at A = 1 and otherwise Eq. (104) shows 
that the summation term must evaluate to (1/2) ln(2/7r). In fact, we may demonstrate 
Lemma 2. For Re A > 0, Re x > 0, and Re y > we have (a) 
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and (b) 



E 



,!-() \ ^ + ^/ ( I) 



a; u r Ut 



v y Jo lnw + 



(A + 1) £(-!)" In 



(A + l) 



m=0 

In — 7 — — r — In 



y + m 
x + m 

r(^) 



(107) 



r(y/2) r(x/2) 
For the proof of part (a), we proceed as in Proposition 8, defining the integral 



(108) 



I(q,z,a) = [ u q 1 (l-au) z l du. 
Jo 



(109) 



We then evaluate 



I(q, z, a)dq 



1 (u y - x 



u 



[1 — au 



,2-1 



du 



(110) 



also by means of binomial expansion of the integrand. Putting z — 1 = n > an 
integer then gives the first part of the Lemma. For part (b) we first use the result 
of part (a), interchanging summation and integration and obtaining Eq. (106). If we 
expand the integrand factor 1/(1 + u) of Eq. (106) as a geometric series, valid for 
\u\ < 1, and then use a tabulated integral [17] (p. 543) we find Eq. (99). Equation 
(108) may be found directly from a known integral [TTJ (pp. 543 or 544) applied to 
Eq. (106), or by using the Hadamard product representation of the Gamma function 
in conjunction with Eq. (107). For the latter we note (cf. [T7], p. 936) 

r r(x/2) 



n i- 



T{y/2) r (2±i) t=o\ (2k + y + l))\ 2k + x 



1 + 



111) 
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Therefore we obtain 

r r(x/2) _ 



In 



5> 



'2£; + y\ (2k+ l + x" 



E 



r( y /2) r (2±i) ^ v^ + ^y v 2A; + 1 + W ,=o 

£(-irin( m+2/ 



In f 2A: + ^ _ ln / /2A; + 1 + ^ N 



2£; + x 



m=0 



2A; + 1 + x. 



(112) 



,m + x , 

and the Lemma is again completed. Alternatively, we may directly relate the left 
side of (106) to the right side of (107) simply by reordering the double summation as 

n=0 2^£=0 — 2^£=0 l^n=i- 

Based upon a very special case of Lemma 2(a) we have 
Corollary 15. We have (a) for Re y > — 1 

"i (u y - 1) 



-du = ln(y + 1), 



(b) 



and (c) 



g(-i)'(*)>G + i) = £ 



J J XJ 

k k 



du 
lnw' 



(113) 



(114) 



s^s^C)^^^^^ 1 ^©- (n5) 



For part (a), we put n — 0, x — 1, and y — > y + 1 in Eq. (105). For part (b), we use 
the integral representation of part (a) and evaluate the two binomial sums. For part 
(c) we find, using part (b), 



A \ k k fk\ 1 00 r°° 

\ E(-1) J • Wjln(j + 1) = £/ TTTnl 



! r ( A _„)»_ (A _r' "'" 



(A+i) /vizjt) *L = i ( A + i) h riv 

2 Jo VI + uj Inu 2 V ; VW 



Inu 
(116) 
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Of many ways to evaluate the last integral of Eq. (116), one may use [T7] (p. 542). 
The Corollary is demonstrated and Eq. (104) is affirmed. 

Moreover, we must recover the values £(— n) = (— l) n B n+ i/(n + 1) for n an integer 
from Eq. (84). For n even this includes the trivial zeros of the zeta function, whereby 
-E>2m+i = for m > 1. We have from Eq. (84) 

«-» - (irW) (TTa, £ (i^) §<-')< + 1)n (117) 

For A = 1 this equation recovers the old formula of Worpitsky for the Bernoulli 
numbers [7J [35] . 

As a byproduct of this work we obtain interesting infinite series (or products) 
for fundamental constants such as the Euler constant and hi7r. The rapidity of 
convergence may make some of these suitable for computational applications. We 
omit many such binomial summations that may be obtained by methods very close 
to Propositions 8 and 9. 

With respect to the right side of Eq. (12), integers m that are powers of 2 play a 
special role. It is only these contributions for which 1 — A(m) = 1 — In 2 > 0, while 
all other integral powers of primes give 1 — A(m) < 0. 

Within the Li criterion for the Riemann hypothesis [26], the sum S^n) is given 
by (e.g., HI]) 

SM = ± (j ^ [ff Ms - DCWU • (us) 

Therefore, by using the Hadamard product representation of the Riemann zeta func- 
tion, it is easy to see that S^n) is connected with sums over its nontrivial zeros. If 
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S^n) has linear or sublinear growth in n, the Riemann hypothesis holds. According 
to our decomposition (15), the validity of the Riemann hypothesis is equivalent to the 
sum S\(n) having linear or sublinear growth in n. In fact, we conjecture (see below) 
that |Sa(^)| = 0(n 1 / 2+t ) for e > arbitrary (This conjecture is stronger than the 
Riemann hypothesis itself.) 

In regard to Proposition 4, we put 

A 2 (x) = d ( n ) - x In x - (2 7 - l)x = 0(x Q2+£ ), e > 0, (119) 

n<x 

wherein o; 2 is the least such number for every positive e. Dirichlet knew that a 2 < 1/2 
and the best result to date may be a 2 < 131/416 [21J. In fact, Huxley showed 
that A 2 (x) = 0(x 23 / 73 In 461 / 146 x) [20] and improved this very recently to A 2 (x) = 
0(x 131/416 ) [2T] . The smallest possible value of a 2 is 1/4, and we prove that \S 7 (n) + 
n\ = 0(n 1//4 ) (see after Eq. (120) below). The method of Propositions 1 and 4 is 
easily extended to many other pole-free combinations of Dirichlet series. 

Proposition 6, its special case Proposition 6.1 of Ref. [9], or other series repre- 
sentations of the Stieltjes constants may be used to obtain alternative summation 
representations of the 0{n) sum S^(n). 

Numerical investigations indicate that S" 7 (n) is close to — n together with a small 
oscillatory component, while S/^(n) is close to n with a small oscillatory component 
5 , 2a(^) [13] . Therefore, the crucial sum 5" 2 (n) appears to arise from substantial can- 
cellation of 0(n), leaving a slowly growing, oscillatory contribution. A demonstration 
that S 2 (n) satisfies a one-sided subexponential bound would suffice to verify the Rie- 
mann hypothesis. 
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As a point of emphasis, the Riemann hypothesis will only fail if a Li/Keiper 
constant becomes exponentially large in magnitude and negative. In particular, 
the Criterion (c) of Ref. [B] now carries over to the crucial subsum S2\(n). Therefore 
the Riemann hypothesis is invalid only if this sum becomes negative and exponentially 
large in magnitude for some n. We may spell this out in the following way. 
Condition. Suppose that there is a value of 1/2 <p< oo such that |<S2a( w )| = 0(n p ). 

Then upon Condition the Riemann hypothesis will follow as a Corollary. It is 
compelling within the Li criterion approach that the optimal order of the sum S2\(Ti) 
is not necessarily required. As indicated, we suspect that the lowest possible order of 
this sum is close to O^n 1 ^ 2 ). 

From Fejer's formula for the asymptotic form of L"(x) we have for n — > oo 

L\_ x {x) = -^e x/2 x- 3/ \n - 1) 1/4 cos(27(n - l)x - 3vr/4) + 0(n- 1/4 ) } x > 0. 
V 71 " 

(120) 

Therefore we now show that the oscillatory component of Sj(n) + n grows as Ofn 1 / 4 ). 
Indeed, the last two terms on the right side of Proposition 2(d) contribute at 0(n 1//4 ), 
with cos(2a/^ — 1?/ + 0) factors. The last three terms on the right side additionally 
contribute at the next lowest order of n _1//4 . For the remaining term on the right side 
of Proposition 2(d) we have 

- J2 L "- l(lnZ/) =-n + J2 L "- l(lnt/) =-n + C(iV 1 /^ n V4 )) (m) 

u=l V u=2 V 

with 5 > 0. We may therefore summarize as 

Corollary 16. As n — >• oo we have \S 7 (n) +n\ — (^(n 1 / 4 ). 
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Besides the indications given in Ref. [TT] that the Laguerre calculus is pervasive 
within the Li/Keiper formulation of the Riemann hypothesis, we have very recently 
systematically presented the structural origins of this framework [12J. The Li/Keiper 
constants arise as a sum over complex zeta zeros of a Laplace transform of the asso- 
ciated Laguerre polynomial L^_ 1 (x). We have 



L n (p) = J™ e-<™-L x n _ x {u)du = 1 - fl - ~) 



(122) 



that vanishes for pj = (1 — e 27nj ) 1 , with j = 1, . . . ,n — 1. These Laplace transform 
zeros have real part 1/2: 



1 (1- e - 2nij / n ) 1 
Pj ~ 2 [1 - cos(2vrj/n)] ~ 2 



1 + i cot 



7TJ 

n 



(123) 



That these quantities lie on the critical line may be more than just a curiosity and 
may partly explain the distinguished role of the polynomial Lj l _ l . 
Furthermore, we recall that [H] (Appendix I) 

L n (s)= f 1 x s - 1 L 1 n _ 1 (-\nx)dx = - 2 F 1 (l-n,l;2;-), (124) 
Jo s V sj 

in terms of the Gauss hypergeometric function. We put C n (s) = s n L n (s) and have 
Corollary 17. We have the functional equation 

L n (s) = -(l-J) L n (l-s), (125) 

and 

C n (s) = (-l) n+1 C n (l-s), (126) 
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that follows immediately. Equation (125) may be obtained directly or by applying 
the transformation formula [T7] 



to the right side of Eq. (124). 

Propositions 1-5 and other results that we have obtained help to expose more 
of the analytic structure of the Stieltjes and rjj constants. We have obtained novel 
integral and other representations of the Stieltjes constants that enable new integral 
and series representations of the sums 5 7 (n) and S2 (n). The growth behavior of 
5*2A(n) and S^n) have direct implication for the validity of the Riemann hypothesis. 




z-1 



z 



) 



(127) 
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